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ABSTRACT 

Accretion and merger shocks in clusters of galaxies are potential accelerators of high-energy protons, which 
can give rise to high-energy neutrinos through pp interactions with the intracluster gas. We discuss the pos- 
sibility that protons from cluster shocks make a significant contribution to the observed cosmic rays in the 
energy range between the second knee at ~ 10 17 ' 5 eV and the ankle at ~ 10 eV. The accompanying cumula- 
tive neutrino background above ~PeV may be detectable by upcoming neutrino telescopes such as IceCube or 
KM3NeT, providing a test of this scenario as well as a probe of cosmic -ray confinement properties in clusters. 
Subject headings: acceleration of particles — cosmic rays — galaxies: clusters — neutrinos 



1. INTRODUCTION 

Clusters of galaxies (CGs) represent the lar gest gravi - 
tationally bound objects in the universe (e.g., IVoit 2005). 
According to standard, hierarchical scenarios of cosmological 
structure formation, they are the latest systems to virialize 
and continue to grow through merging and accretion of dark 
matter and baryonic gas, thereby generating powerful shock 
waves on Mpc scales. In particular, accretion shocks with 
high Mach numbers are expected on the outskirts of massive 
CGs, potentially leading to efficient acc e leration of high - 
energy particles (e.g., iMiniati et al. 2000t iRyu et al. 20031) . 
Here "accretion" signifies not only infall of diffuse intergalac- 
tic gas, but also minor merger events that induce sufficiently 
strong shocks near the virial radii. Moderate Mach number 
shocks arising further inside t he CG could al so be important 
in certain situations (e.g., iRvu et al. 2003b . Predictions 
for the associated nonthermal radiation, notably high- 
energy gamma-r ays, have been d isc ussed by a number of 
authors (e.g.. IVolk et al. 19961 IBerezinskv et al. 1997|: 
Colafrancesco & Blasi 1998; lLoeb & Waxman 2 00C 



Inoue et al. 2005tlAndo & Nagai 2 008 ) 



Cosmic rays (CRs) are observed over 1 1 decades of energy 
from ~ 10 9 eV to > 10 20 eV, and their origin is under intense 
debate. The all-particle spectrum is characterized by broken 
power-laws with a number of breaks: the knee at ~ 10 15 5 
eV where the spectral index changes from ~ -2.7 to ~ -3.0, 
the second knee at ~ 10 17 5 eV where it changes from ~ -3.0 
to -3.2, and the a nkle at ~ 10 18 5 eV wher e it changes from 
~ -3.2 to ~ -2.7 (iNaeano & Watson 200(1 . Galactic super- 
nova remnants (SNRs) are widely believed to be responsible 
for CRs at least up to the knee, a nd probably up to some- 
what higher energies dHillas 2 005. and references therein). 
In contrast, ultra-high-energy CRs (UHECRs) with ener- 
gies abov e the ankle are generally t hought to be extragalac- 
tic (e.g.. lNagano& Watson 20001: IGaisser & Stanev 20061 
Bergman & Belz 2007). The Pierre Auger Observatory has 
revealed an anisotropy in the arrival directions of UHE- 
CRs that supports this view (Abraham et al. 2007; 2008; 
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see however Abbasi et al. 2008) but the actual iden- 
tity of the sources remain unknown (e.g., iKachelrieB 20081: 
llnoue 20081) . The origin of CRs with energies between the 
second knee and the ankle is even more controversial, and 
both extragalactic and Galactic viewpoints have been advo- 
cated. _TJie_sourc^s_jnaybe th£_sam£_as_those for UHE- 
CRs (IBerezinskv et al. 20061: lAloisio et al. 2007b . or some- 
thing entirely different such as speci fic types of supernovae , 



hypernovae, Galactic winds, etc (IGaisser & St anev 2006 
iBergman & Belz 20071; IWang et al. 20071: iBudnik et al. 2008 , 
and references there in). 

In this letter, we discuss high-energy neutrino emission 
from CGs in light of the possibility that accretion and/or 
merger shocks in CGs are important sources of CRs with 
energies between the second knee and the ankle. The 
observability by n ear future tele scopes suc h as IceCube 
dAhrens et al. 20041) or KM3NeT dKatz 20061) is addressed. 
Note that this view is different from scenarios in which clus- 
ter shocks are the main sources of UHECRs above the an- 
kle dNorman et al. 19951; iKang et al. 19961; iKang et al. 19971; 
llnoue et al. 20071) . More detailed discussions including the 
accompanying hadronic gamma-ray emission will be pre- 
sented in a subsequent paper. We adopt the ACDM cosmo- 
logical parameters ft/, = 0.04, fl m = 0.3, fl\ = 0.7 and h = 0.7. 

2. COSMIC-RAY PRODUCTION 

We first estimate the maximum energy of the acceler- 
ated CRs. The virial radius of a CG with mass M = 
10 15 M 15 M Q is r vir ~2.4MpcM 1 1 5 /3 J F(z,^ m )(/!/0.7)" I (l+z)- 1 , 
where F(z, Sl m ) is a fact or of order unity that depends weakly 
on redshift z and Sl m (IVoit 200 5). We write the shock ra- 
dius as r s h = As h^vir, with A s h ~ 1-10 expected for accre- 
tion shocks (e.g jRyu et al. 20031) . The typical shock velocity 
V s h should be comparable to the velocity of the infalling gas 

Vff ~ 2000 kms" 1 Mj 5 /2 (r sh /l Mpc)" 1/2 (llnoue et al. 2005b . 

The typical shock acceleration time for CRs with energy 
e and cha rge Z is t» cc « 20K, h M = (20/3)(ce/ZefiV 2 ) £ 
(see, e.g., iBlandford & Eichler 19871 for reviews). Here B 
and n s h are respectively the magnetic field and diffusion co- 
efficient at the shock, and £ = (B/ SB) 2 where £ — > 1 in the 
Bohm limit. Although the magnetic fields at cluster shocks 
are uncertain, we take B ~ 1 pG, as supported by recent 
X-ray ob servations of diffuse radio relics near r v j r for sev- 
eral CGs dFerreti & Neumann 20061 iHubert Chen et al. 20081; 
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iNakazawa et al. 20081) . We also post ulate £ ~ 1, as ob - 
served to be the case for some SNRs (Uchiy amaetal. 20"07h . 
and which can be expected if the fields are generated lo- 
cally at the shocks by mec hanisms such as the CR stream- 
ing instability (see, e.g.. iBeil & Lucek 20011: iBell 2004t 
Vladi morov et al. 20 06: Bvko v et al. 2008b . 

The maximum energy of the accelerated CRs £ max can 
be estimated by equating f acc with various limiting time 
scales, such as the diffusive escape time from the shock 
fesc ~ ff h /6K s h, and the energy loss time due to photo- 
hadronic and/or photodisintegration inter actions with the 
CMB and the in frared (IR) background dKangetal. 19971: 
llnoue et al. 20071) . When the shock is due to a transient 
merger-like event, the lifetime of the shock may also be rel- 
evant, which can be estimated by the dynamical time fdyn = 
'"sh/Vsh- In me latter case, e max would generally be deter- 
mined by f dyn so that e max ~ 1.6 x lO^eVZ^Mjf S-eA^ 2 , 
where B = B-^pG. Hence, we expect that cluster shocks 
can accelerate at l east protons up to around the ankle 
dNorman et al. 1995b . 

Next we consider the energetics. Assuming a total mass 
accretion rate M rj O.IV^/G and gas fraction f g = f2/,/T2 m ~ 
0.13, the dissipation rate of infalling gas kinetic energy 
through the accretion shock of a CG with mass M is estimated 

to be L ac » f e GM M/r sh ~ 7 x lO^ergs" 1 (/ ? /0.13)Mi 5 /3 
dKeshet et al. 20041) . Taking the local density of massive 
CGs n CG (M > Mis) ~2x 10- 6 Mpc" 3 dJenkins etal. 2001b 
and a CR injection efficiency e acc ~ 0.2, the CR power 
from CGs per logarithmic energy interval at e = 10 I8 eV 
is e 2 (dri/de) ~ lO^ergMpc^yr" 1 (R/50)' 1 . Here R(e) = 
(f^de'e'(dN/de')j /(e 2 dN/de) depends on the injection 

CR spectrum; in the case of a single power-law with index p 
and minimum energ y = 1 GeV, R ~ 25 for p = 2.0, and 
R ~ 300 for p = 2.2 (favlurase et al. 2008b . 

In comparison, the observed CR spec- 
trum for 10 17 eV < e < 10 18 ' 5 eV is $ ~ 
9.23 x 10" 28 ey ~ 1 cm- 2 s- 1 sr 1 (£/6.3 x 10 18 eV)" 3 ' 2 
(Nagano & Watson 2000). The implied CR source power 
and spectrum depend on evolution effects and the uncertain 
CR composition at these energies . Here we assume that 
they are proton-domina nt as in some extragalactic scenarios 
dBerezinskv et al. 2006b . and cases with a more general com- 
position will be discussed in the future work. The required CR 
power at e = 10 18 eV is then s 2 (dri/de) ~ 10 45 ergMpc~ 3 yr _1 , 
crudely accounting fo r energy losses during diffusive inter- 
galactic propagation dBerezinskv et al. 2006b . As long as 
R < 100, the two powers would be comparable and CGs 
energetically viable. However, since intergalactic propa- 
gation should steepen the spectral index fro m the injection 
value p by ~ 0.5-0.8 dBerezinskv et al. 2006b . the observed 
index of -3.2 requires p ~ 2.4-2.7 and hence R ;§> 100 for 
single power-law spectra extending down to £ m ; n = 1 GeV. 
Avoiding excessive energy demands motivates a broken 
power-law form with break energy £/,: dN/de oc e~ p ' for 
£ < Eb and dN/de oc e~ P2 e\p(-(e / e max )) for e > e b (see also 
lAloisio et al. 2007b . 

We consider two possibilities as to how such spec- 
tra may actually occur. One is through the superposi- 
tion of hard spectra (pi ~ 2 .0) with a distribution of £ max 
dKachelrieB & Semikoz 20 06). which can be related to ac- 
cretion shocks with a distribution of M. It was seen that 
£max oc M 2 / 3 if the relevant condition is f acc w fd yn and if the M- 



dependence of B is not strong. A realistic CG mass function 
«cg(> M) ocM _1 exp(-(M/1.8 x 10 14 M Q )) can be approxi- 
mated over a limited range of M as a power-law «cg(> M) oc 
M~ a dJenkins et al. 2001b so that dN/de oc £~^~5 for e > e b . 
This could allow p\ ~ 2.0 and p2 ~ 2.0-3.3. 

Another possibility is a two-step acceleration process, a 
first source providing a seed CR population with hard spec- 
tra (p\ ~ 2.0) up to eb, which is then picked up by a sec- 
ond source and accelerated further with softer spectra to £ max . 
Since CRs with sufficiently low energies are likely to be con- 
fined in the intra c luster medium (ICM) f or very long times 
dVolketal. 19961: iBerezinsky et al. 1997b . the seed popula- 
tion can come from a number of sources, all accumulated 
over the history of the CG: the low energy portion of accre- 
tion shock CRs, supernova-driven galactic winds (GWs), and 
the jets of radio-loud active galactic nuclei (AGNs). Their 
relative importance can be estimated through their contribu- 
tions to the heating of the ICM, which should be roughly pro- 
portional to their CR output as long as the relevant shocks 
are sufficiently strong. In the absence of GWs or AGNs, 
high Mach number accretion shocks are expected to con- 
tribute ~10% of the heating of the ICM, while the re- 
maind er is mediated by low Mach number merger shocks 
(e.g., iRyu et al. 2003b . GWs are unlikely to play a sig- 
nificant role in ICM heating due to severe radiativ e losses 
during their formation (e.g. Krav tsov & Yepes 200 0). (Note 
that CRs of supernova origin escaping from within the host 
galaxy will also suffer heavy adiabatic losses). In con- 
trast, AGN jets can contribut e 1-2 keV per bary on of heat 
input directly to t he ICM (|ln oue & Sasaki 2001) (see also 
lEnBlin et al. 19971: lEnBlinetal. 1998b . For massive clus- 
ters with temperatures ~10 keV, this implies that CRs from 
AGNs can be energetically comparable to those from accre- 
tion shocks, and may be even higher for less massive clus- 
ters. Subsequent acceleration of these seed CRs to £ max 
may be achieved through merger and/or accretion shocks 
with moderate Mach numbers M. ~ 2.5-5, leading to p2 ~ 
2.2-2.7. The break energy e b may correspond to the con- 
finement energy e^s above which CRs begin to escape dif- 
fusively out of the ICM. Under Kolmogorov-like turbulence, 

e m - 1.8 x 10 17 eVZ(r/1.5Mpc) 6 Kc 3 G! 3o( Af / 1 Gyr)~ 3 , where 
kqg = 10 30 kcg.30 cm 2 s _1 is the diffusion coefficient in the 
ICM at £ = 1 GeV, and At is the time el apsed after injection 
dVolket al. 19961: IBerezinsky et al. 1997h . Within the uncer- 
tainties, we see that CGs could be a viable source of CRs with 
energies between the second knee and the ankle. 

3. NEUTRINO PRODUCTION 

We now evaluate the spectra of associated gamma-rays 
and neutrinos, which are inevitably gene rated through pp 
interactions with the a mbient ICM gas dVolketal. 1996t 
IBerezinsky et al. 1997t IColafrancesco & Blasi 19981) . In 
view of the above, we assume that CRs with a broken 
power-law spectrum is realized with p\ = 2.0 and p2 = 
2.4. We choose e b = 10 165 eV or 10 175 eV, giving re- 
spectively R ~ 78 or 35. The spatial distribution of the 
thermal ICM g as is generally well-con strained from X-ray 
observations dPfrommer & EnBlin 20041) . However, that for 
the CRs is uncertain, and we consider the following four 
models. Model A: CRs are uniformly distributed within 
r s h, with A s h chosen such that fd yn = 1 Gyr. Model B: 
CRs are uniformly distributed within r v ; r , giving a conser- 
vative estimate compared to other models. Isobaric: CRs 
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Fig. 1. — Expected event rates for muon neutrinos {Vp + Pp) in 
IceCube-like detectors from five nearby CGs: Virgo, Centaurus, 
Perseus, Coma and Ophiuchus. Broken power-law CR spectra with 
pi = 2.0, f2 = 2.4 and ej, = 10 17 5 eV is assumed, and the Isobaric 
model with X C r = 0.029 is used. Note that IceCube and KM3NeT 
mainly cover the northern and southern celestial hemispheres, re- 
spectively. Neutrino oscillation is taken into account. 
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Fig. 2. — Cumulative neutrino (v e + v e + Vp, + Dp + v T + v T ) back- 
ground from CGs for broken power-law CR spectra with pi = 2.0 
and p2 = 2.4. The break energies are ej = 10 17 ~ eV (thick lines) and 
e b = 10 16 5 eV (thin lines), respectively. The CR power is normal- 
ized to e 2 (dri/de) = 2 x 10 45 ergMpc~ 3 yr~' at e = 10 18 eV, as required 
to account for CRs above the second knee. For isobaric, the corre- 
sponding Xqr is 0.029 and 0.067. For central AGN, Kolmogorov-like 
turbulence is assumed with kcg = 10 30 cm 2 s~' . We take td yn = At = 1 
Gyr and z max = 2. WB represents the Waxman-Bahcall bounds 
dWaxman & Bahcall 19981) . 

at each radii have energy densit y proportional to that of 
the thermal gas wit h ratio Xqr ( Pfrommer & Enfili n~2004t 
lAndo & Nagai 2 008). Central AGN: CRs are distributed 

as dN/de oc r _I £"''"5 for e > (r 2 /6«;cGAf) 3 and dN/de oc 

e~ p for e < (r^/6KcGAf) 3 , corresponding to CRs diffusing 
out from a central source such as an AGN as discu ssed in 
(iBerezinsky et al. 1997t IColafrancesco & Blasi 19981) : here 
accretion or merger shocks may not be involved. But, in the 
two-step acceleration scenario, this model could be more re- 
alistic below the break energy. We perform numerical calcu- 
lations of the neutrino spectra using formulae ba sed on the 
SIBYLL code at high energies dKelner et al. 20061) . 
The neutrino and gamma-ray fluxes can be estimated an- 



alytically via the effective optical depth for the pp reaction 
as f pp m 0.8tT w ,«NCfint, where «n is the target nucleon den- 
sity in the ICM, a pp is the pp cross section and t mt ~ fd yn 
or max[r/c,/diff] is the pp interaction time. Because «n ~ 
lQ-^cnT 3 at r ~ 1.5 M pc (IColafrancesco & Blasi 19981: 
iPfrommer & Enfilin 2004 and a pp ~ 10" 25 cm 2 in the 100 
PeV range dKelner et al. 20061) . we obtain 

fpp -2.4 x 10- 3 n N ,-4. 5 (Wl Gyr). (1) 

Roughly speaking, high-energy neutrinos from charged-pion 
decay have typical energy e„ — 0.03e (true only in the aver- 
age sense, because charged particles have wide energy dis- 
tributions an d high multipliciti es as expected from the KNO 
scaling law) dKelner et al. 20061) . Hence, neutrinos > PeV are 
directly related to CRs above the second knee. 

First we obtain numerically the neutrino spectra and ex- 
pected event rates from five nearby CGs, utilizing the model 
or double-/? model description in Tables. 1 and 2 in Pfrommer 
& EnBlin 2004 for the thermal gas profile of each CG (Fig. 
1). Our gamma-ray fluxes for single power-law spectra agree 
with the results of Pfrommer & EnBlin 2004. As is apparent 
in Fig. 1, the detection of neutrino signals from individual 
CGs could be challenging even for nearby objects. It may be 
achievable, however, through a detailed stacking analysis. 

More promising would be the cumulative background sig- 
nal. A rough estimate of the neutrino background is (e.g., 
iMurase 20071: IWaxman & Bahcall 19981) 
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where CGs are assumed to be the main sources of CRs from 
the second knee to the ankle. Here, «cg(0) is the local density 
of mas sive CGs and f, is a correction factor for th e source evo- 
lution dMurase 2 007: Wa xman & Bahcall 19981) . For detailed 
numerical calculatio ns of the background, we trea t more dis- 
tant CGs following (Colafrancesco & Blasi 1998), adopting 
the mass function of d Jenkins et al. 200lh . The results for 
the broken power-law case are shown in Fig. 2. With et, = 
10 17 5 eV, the expec ted event rates above 0. 1 PeV in IceCube 
dAhrens et al. 20041) are ~ 2 yr" 1 for model A, — 1 yr" 1 for 
model B, ~ 5 yr -1 for isobaric and — 3 yr" 1 for central AGN. 

Hence, upcoming telescopes may be able to find multi- 
PeV neutrino signals from CGs, providing a crucial test of 
our scenario. From Eq. (0, we can also estimate the cor- 
responding gamma-ray background from ir° decay, which is 
el f <^ 1 — (10 -9 - 10" 8 )GeVcm" 2 s" 1 str" 1 for the broken power- 
law case. This is only (0.1 - 1)% of the EGRET limit, con- 
sistent with the nondetection so far for individual CGs. Note 
that the expected gamma-ray background flux would increase 
if Eb can be decreased, requiring larger CR power from CGs. 

4. IMPLICATIONS AND DISCUSSION 

To test the CG origin of second knee CRs, high-energy neu- 
trinos should offer one of the most crucial multi-messenger 
signals. Unlike at the highest energies, CRs themselves in 
the 10 18 eV range offer no chance of source identification as 
they should be severely deflected by Galactic and extragalac- 
tic magnetic fields. Moreover, due to magnetic horizon ef- 
fects, extragalactic CRs < 10 17 eV may not reach us at all 
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dLemoine 20051 : iKotera & Lem oine 2007]), so even the bro- 
ken power-law spectral form will not be directly observable. 
Gamma-rays are unaffected by intervening magnetic fields, 
but those at > PeV energies relevant for the second knee are 
significantly attenuated by pair-creation processes with th e 
CMB and cosmic IR backgrounds (e.g., iKachelrieB 20081) . 
In constrast, neutrinos in the Pe V-EeV energy range should 
be unscathed during propagation dBhattachariee & Sigl 20001 
and references there in). Consequently, such neutrinos may 
also constitute a unique tool for probing the uncertain CR con- 
finement properties of CGs through the dependence on E&g . 

AGNs can complicate the cluster shock neutrino signal, ei- 
ther by emitting PeV-EeV neutrinos themselves, or inject- 
ing CRs that produces neutrinos without the intervention of 
cluster shocks. In principle, cross correlation of the detected 
events with known CGs and AGNs should be an effective dis- 
criminant. In the former case, AGNs inside as well as outside 
CGs should correlate, whereas in the latter, CGs with power- 
ful AGNs should correlate stronger than those without. For 
cluster shock neutrinos, correlation with all sufficiently mas- 
sive CGs is expected. 

Gamma-ray observations at GeV-TeV energies would also 
be crucial. In combination with > PeV neutrino observa- 
tions, they can probe the CR spectrum over a broad en- 
ergy range and test our broken power-law assumption. By 
providing information on the spatial distribution of sub-PeV 
CRs, they would also help to distinguish among our differ- 
ent models, and to constrain the AGN contribution of CRs 
and neutrinos. However, if £j is sufficiently high, their detec- 
tion may not be trivial except for a few nearby CGs such as 
Virgo. Other emission processes may al so be at work (e.g., 
Loeb & Waxman 2000; Inoue et al. 2005), complicating the 
extraction of the tt decay component. More details on the 
gamma-ray emission will be given in a following paper. 



Note that high-energy neutrinos can also be produced 
by photomeson in teractions with the IR background (e.g., 
iTakami et al. 2007|). which may possibly be enhanced inside 
CGs dDe Marco et al. 20061) . The effective optical depth for 
this process is roughly f pl ~ 10~ 3 - 10 _1 . These neutrinos be- 
come important above ~ 0. 1 EeV, where detection by IceCube 
or KM3NeT is relatively more difficult. We defer the study of 
such photomeson-induced neutrinos in CGs to the future. 

Our neutrino predictions were based on simplified models 
for the CR distribution in CGs. More realistic evaluations 
need to include the inhomogeneous structure of accretion and 
merger shocks, the magnetic field distribution inside and out- 
side CGs, and the associated CR acceleration and propagation 
processes. The CR confinement properties are especially cru- 
cial. Although somewhat extreme, if CR confinement in CGs 
is more efficient than we have assumed, the contribution to 
the CR flux would decrease while the relative neutrino flux 
would increase. Confinement may possibly be effective even 
on larger scales such as filaments or superclusters, in which 
case the CR contribution could increase relative to the neu- 
trino flux, as our Galaxy resides in the local supercluster. 
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